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1. Introduction 

In November 1915, Albert Einstein published the final paper that eompleted his theory of 
gravity and spaeetime now known as general relativity (GR), ehanging our view of gravity and 
spaeetime for ever. During its first fifty years, GR was eonsidered a theoretieal tour-de-foree but 
one with meager observational evidenee. Only in the 1960s did teehnology usher in the remarkable 
field of experimenfal gravity (Misner el al. 1973; Will 1993). Over Ihe subsequenf fifty years, GR 
has passed all experimenfal fesfs in laboralories, in Ihe Solar Syslem and in various slellar sys¬ 
tems, parlieularly in binary pulsars (Will 2014). Aslrophysieal observalions using radio bands have 
played one of Ihe mosl imporlanl parls in Ihe history of lesling GR, wilh Ihe abilily to make preeise 
measuremenls and probe info Ihe vieinily of eompael bodies, namely blaek holes (BHs) and neu- 
Iron slars (NSs). Parlieularly, radio liming of binary pulsars has preeisely probed Ihe gravifalional 
properlies of NSs and, for Ihe firsl lime, lesled Ihe radialive properties of gravily, demonslraling 
lhal gravifalional waves exisl and lhal eompael binaries lose energy from Iheir emission al Ihe rales 
predieted by GR (Hulse & Taylor 1975; Taylor el al. 1979; Kramer el al. 2006; Weisberg el al. 
2010; Freire el al. 2012b; Anloniadis el al. 2013). This eonfirmalion and preeise eharaelerisalion 
of Ihe radiative properties of gravity is nol only erilieal from Ihe poinl of view of our underslanding 
of fundamenlal physies and many aspeels of aslrophysies, bul il will also open, via gravifalional 
wave deleelors, a whole new window on Ihe Universe. The Square Kilomelre Array (SKA), as Ihe 
mosl powerful member of Ihe nexl generation of radio leleseopes, will allow mueh more preeise 
radio timing of pulsars and will diseover dozens of rare relalivislie binary systems. For Ihis reason 
il will play a unique role in many areas of experimenfal gravity. 

Despite ils sueeess over Ihe Iasi eenlury, GR still faees greal ehallenges on Ihe researeh fron¬ 
tiers of modern physies. Firsl, il is known to be ineomplele, sinee il fails al Ihe eenlre of BHs. 
Seeond, as a non-renormalisable field Iheory, il appears to be ineompalible wilh quanlum prinei- 
ples. Reeoneiling gravity wilh quanlum meehanies slands oul as one of Ihe greal ehallenges of 
fundamenlal physies today. Furlhermore, eoneepls like dark mailer and dark energy pose addi¬ 
tional ehallenges for our underslanding of gravity, and modifiealions of gravity Iheories have been 
proposed to explain Ihese phenomena (Capozziello & de Faurenlis 2011), like Ihe Tensor-Veelor- 
Sealar (TeVeS) Iheory, whieh was advaneed to explain dark mailer (Bekenslein 2004). 

The imporlanl Ihing to us is lhal pulsar timing ean be used to lesl some of Ihese Iheories. Pulsar 
timing uses large radio leleseopes to reeord Ihe times of arrival (TOAs) of pulsed signals, produeed 
by Ihe relation of Ihe pulsar. These TOAs depend on Ihe rolalional and aslromelrie parameters of 
Ihe pulsar, dispersion in Ihe interstellar medium lhal Ihe signals Iraverse, and Ihe motion of Ihe 
radio leleseope in Ihe Solar System. If Ihe pulsar is in a binary, Ihe TOAs also depend on Ihe orbilal 
dynamies of Ihe binary, whieh are determined by Ihe underlying gravifalional Iheory (Weisberg 
& Taylor 1981; Damour & Taylor 1992; Edwards el al. 2006). Some Iheories of gravity, like a 
family of TeVeS-like Iheories inlrodueed in Freire el al. (2012b), resull in orbilal dynamies lhal 
is suffieienlly differenl from GR to be deleeled in Ihe timing of some binary pulsars. Therefore, 
preeise timing measuremenls of Ihese binary pulsars ean eonslrain, or allogelher rule oul, some of 
Ihese gravity Iheories. 

For example, in Ihe Jordan-Fierz-Brans-Dieke Iheory (Jordan 1959; Fierz 1956; Brans & Dieke 
1961), gravity is mediated by a sealar field, (p, in addition to Ihe eanonieal melrie field, g^y. The 
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Figure 1: Left: Parameter space for quantifying the strength of a gravitational field (Psaltis 2008). The 
horizontal axis measures the gravitational potential and the vertical axis measures the spacetime curvature of 
an orbit with a semimajor axis, a, and a total mass, M (here G is the gravitational constant, and c is the light 
speed). Solid and empty squares are known neutron star - neutron star (NS-NS) and neutron star - white 
dwarf (NS-WD) binaries, respectively. The diamond is a known neutron star - main sequence star binary. 
Circles from large to small are pulsar - Sgr A* BH, pulsar - 5000 Mq BH, and pulsar - \QMq BH binaries 
with their assumed orbital periods aside. Right: The compactness of the gravitating companion star versus 
the compactness of the orbit; Me and Rc are the mass and the radius of the companion star, respectively 
(Kramer et al. 2004). 


variation of q) in spacetime introduces observable effects in gravitational experiments, like the time- 
variation of the local gravitational constant G and the existence of gravitational dipole radiation in 
a binary. In a more general class of scalar-tensor theories (Damour & Esposito-Farese 1993, 1996), 
the coupling strength of the scalar field (p to matter depends on the field ifself ip as ao + <Pj8o, where 
«() and jSo fully characterise fhe fheory. Wifhin some of fhe (aojjSo) space, non-perfurbafive effeefs, 
called “sponfaneous scalarisafion”, may happen in compacf bodies like NSs (Damour & Esposifo- 
Farese 1993), which can influence fhe orbifal dynamics of a binary system dramatically. Similar 
effecls are also possible wifhin an exfended family of TeVeS-like fheories (Freire el al. 2012b). Wifh 
fheir abilily fo probe sfrong-field effecls, binary pulsars are ideal laborafories wifh which lo sludy 
Ihese alternative gravilafional fheories. Wifh fhe non-delecfion of gravifalional dipole radiation 
from radio pulsar timing, binary pulsars have already severely conslrained fhe parameter spaces of 
Ihese fheories and provided fhe mosl slringenl lesls of fhe sfrong-field* and radiative properties of 
gravity (Freire el al. 2012b; Anloniadis el al. 2013; Slairs 2003; Wex 2014). 


'Compared to the coalescence of two compact objects (that will be the main source for ground-based gravitational 
wave detectors), two components of a binary pulsar are well separated, that might be prejudicially termed as “weak 
field”. However, as shown by explicit calculations, “strong-field effects” associated with strongly self-gravitating NSs 
could have significant effects on the binary orbital dynamics (Damour & Esposito-Farese 1993). 
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1.1 The role of the Square Kilometre Array 

The Square Kilometre Array will have a major impaet on radio astronomy. The SKA Phase 1 
(SKAl) will feature roughly half of its total eolleeting area within a dense eore, in the eases of 
SKAl-LOW and SKAl-MI D . By synthesising a eoherent sum of these elements, it will be possible 
to deeply seareh for new relativistie binary pulsars, over the whole sky, and to subsequently time 
them with great preeision using the full SKAl-MID array {i.e. ineluding the longer baselines). 
The SKA Phase 2 (SKA2) will provide another major leap in sensitivity, and thus also in timing 
preeision. Both phases of the SKA will signifieantly advanee experimental tests of gravitational 
theories to a new level (Kramer et al. 2004). 

This will happen through two avenues. On the one hand, SKAl-MID and SKA2 will strongly 
improve the quality of the observations of the most interesting binaries eurrently known (e.g., 
PSR J0737—3039 (Burgay et al. 2003; Lyne et al. 2004; Kramer et al. 2006)) and triple systems 
(e.g., PSR J0337-I-1715 (Ransom et al. 2014)), in turn dramatieally improving the quality of the 
existing tests and enabling proposed new tests. Compared with the eurrent status quo, SKAl-MID 
will improve the timing preeision for southern pulsars by about an order of magnitude, and SKA2 
will improve it by up to two orders of magnitude. Different pulsars will gain different improve¬ 
ments in preeision, depending on their sky loeations, pulse eharaeteristies, timing stabilities, and 
so on. 

For preeision pulsar timing, raw sensitivity is one of the most important faetors. High eadenee 
may also be needed for some experiments, in whieh ease sub-arraying will be needed. Typieally, 
TOAs for a partieular pulsar are eolleeted every few weeks, though in eome eases a mueh higher 
eadenee is neeessary, e.g. to foeus on orbital eonjunetion in order to measure the Shapiro delay. 
Sub-arraying ean also help if many sourees need to be timed, and if redueed sensitivity is aeeept- 
able. Though SKAl-MID will provide the highest-preeision TOAs, regular measurements with 
SKAl-LOW may also be needed to traek the time-variable interstellar propagation effeets that also 
influenee the TOAs. In all eases, eoherently dedispersing the data is essential for aehieving the 
maximum timing preeision and aeeuraey. Sueh timing observations ean already begin in an early 
phase of SKAl (say, with 50% of the eolleeting area available), though they will naturally obtain 
less preeise TOAs than with the eomplete array. For some sourees 50% sensitivity will be suffi- 
eient, while for others eertain relativistie effeets may only beeome elearly deteetable onee the full 
array is in plaee. 

The SKA will also discover new relativistie systems with whieh to earry out similar and eom- 
pletely novel tests. With a larger field of view and high sensitivity, the SKA is more effieient in 
finding new pulsars. The larger sensifivify of fhe SKA means fhaf a shorfer infegrafion lime is 
needed lo aehieve fhe required signal lo noise ratio. Consequenlly, in an “aeeeleralion seareh” wilh 
an assumed eonslanl orbilal aeeeleralion (Lorimer & Kramer 2004), fhe smearing due lo varying 
aeeeleralion in an aelual binary will be redueed, making Ihe SKA mueh more effeelive lhan eurrenl 
leleseopes al finding relativistie binaries. Simulations show lhal Ihe SKAl ean diseover a lolal 
of aboul 10,000 normal pulsars and perhaps as many as 1,800 milliseeond pulsars (MSPs), wilh 
SKAl-LOW surveying Ihe sky wilh Ihe Galaelie lalilude \b\ > 5°, and SKAl-MID surveying Ihe 
sky wilh Ihe Galaelie lalilude \b\ < 10° (Keane el al. 2014). The SKA2 will diseover even more 
pulsars, and if is possible lhal among Ihem Ihere will be pulsar - blaek hole (PSR-BH) binaries 
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(Kramer et al. 2004; Keane et al. 2014). The discovery of even a single tight PSR-BH system will 
open a new era of studying BH physics with great precision, including possible tests of the “cos¬ 
mic censorship conjecture” and the “no-hair theorem” (Damour & Esposito-Farese 1998; Wex & 
Kopeikin 1999; Kramer et al. 2004; Liu 2012; Liu et al. 2012). Figure 1 shows the vast unprobed 
regions in the gravity sector that can be probed with different kinds of PSR-BH systems. 

2. Relativistic Binaries 

Pulsar timing takes advantage of the tremendous rotational stability of pulsars, allowing us to 
treat them as free-falling clocks. As a pulsar is monitored over months and years, an ephemeris 
is calculated which accounts for every rotation of the pulsar and can be used to predict future 
pulse TOAs. A pulsar ephemeris always incorporates the spin frequency and generally at least one 
spin derivative accounting for rotational energy loss, plus astrometric parameters. Binary pulsar 
orbits require a further translation of reference frame, accomplished for most binary systems by 
the five Keplerian parameters describing an eccentric Keplerian orbit. In order to account for 
relativistic effects beyond a Keplerian orbit, a set of theory-independent “post-Keplerian” (PK) 
parameters are introduced (Damour & Deruelle 1986; Damour & Taylor 1992) (see Table 1 for an 
incomplete collection of the most important PK parameters). In practice, PK parameters could be 
contaminated by astrophysical effects other than relativistic effects due to gravity. For example, the 
observed orbital decay parameter, has “kinematic” contributions, which result from the proper 
motion of the system and the difference of Galactic accelerations between that system and the Solar 
System Barycentre. Here unless otherwise stated, we assume that non-gravitational contributions 
are properly corrected or demonstrably negligible in effect. If so, in the case of two well-separated 
masses with negligible spin contributions, the PK parameters are functions of the well-measured 
Keplerian parameters, the component masses, the equation of state (EOS) of stellar matters, and 
the parameters describing the gravitational theory (Damour & Taylor 1992) (in GR, the internal 
constitution of the star is iiTelevant to a high “post-Newtonian” (PN) order, only the total masses 
matter — this is the “effacement” property to be discussed in section 3). 

Measuring two PK parameters, we can use a specific theory of gravity to determine the com¬ 
ponent masses of the system; these are important to study stellar evolution theories and, in some 
cases, to constrain EOS (Watts et al. 2014). If more PK parameters are measured, then the theory 
can be tested by a self-consistent argument — using the masses derived in the first stage one should 
be able to predict the subsequent PK measurements. Another way of viewing this is the so-called 
“mass-mass diagram” (see Figure 2 for a GR-based example for PSR J0737—3039; Kramer et ah, 
in prep.). For a gravitational theory to pass the test(s), all curves in the diagram should intersect 
in some region, i.e., the theory must be able to describe the component masses in a self-consistent 
way (Taylor & Weisberg 1982). 

The current state-of-the-art binary in this area is the Double Pulsar (PSR J0737—3039A/B) 
(Burgay et al. 2003; Lyne et al. 2004), which enables five fesfs of GR. Figure 2 shows ifs mass- 
mass diagram, where seven curves corresponding fo seven mass consfrainfs infersecf af one poinf 
(Kr'amer ef al. 2006; Breton el al. 2008; Kramer & Wex 2009), indicating lhal GR provides a self- 
consisfenl description of all of ifs liming measuremenfs. The Double Pulsar resides in Ihe soulhern 
sky and will be an imporlanl largel for Ihe SKA. The timing precision of Ihe SKA will enable fesfs 
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Table 1: The most important PK parameters that could be used in pulsar timing of binaries (Damour & 
Deruelle 1986; Damour & Taylor 1992; Wex & Kopeikin 1999; Lorimer & Kramer 2004; Edwards et al. 
2006). In practice, for a specific binary pulsar, only some PK parameters are measured, depending on the 
characteristics of the pulsar timing experiment. 

Parameter 

m 

time derivative of the longitude of periastron ft) 

7 

amplitude of the Einstein delay 

Pb 

time derivative of the orbital period 

r 

range of the Shapiro delay 

s 

shape of the Shapiro delay 

^SO 

precession rate of the pulsar spin 


mismatch in eccentricities (see text) 

e 

time derivative of the orbital eccentricity e 

X 

time derivative of the projected semimajor axis of the pulsar orbit x 

(b 

second time derivative of the longitude of periastron 

X 

second time derivative of the projected semimajor axis of the pulsar orbit 


of relativistic effects beyond their first-order PN approximations (Kramer & Wex 2009). Measuring 
higher-order PN effects is important since it allows the investigation of the effects from pulsar spin 
(Barker & O’Connell 1975; Damour & Schaefer 1988; Kramer & Wex 2009). The fractional 
effects due to the next PN order^ on the advance rate of periastron, namely and on 

the orbital decay, namely p3-5PN^p2.5PN^ order of Simulations showed 

that by obtaining a timing variance below ~ 5ps for the Double Pulsar we will reach the precision 
required to detect higher-order effects (Kramer & Wex 2009); this is very likely with SKAl-MID, 
scaling from current observations. The SKA2 will have sub-/rs timing precision for the Double 
Pulsar (assuming its pulse profile shape can sfill be well modelled), raising fhe prospecf of highly 
precise measuremenfs of fhese higher-order parameters. 

Besides fhe orbifal dynamics of fwo poinf masses, spin confribufions may play an imporfanf 
role in fhe observafions of pulsar binaries (Barker & O’Connell 1975). Due fo fhe curvafure of 
fhe spacefime produced by fhe companion sfar, fhe rofafion axis of a freely falling objecf (here fhe 
pulsar) suffers a precession wifh respecf fo a disfanf observer; fhis effecf is known as “geodetic pre¬ 
cession”. Because of if, over lime differenf emitting regions of fhe magnelosphere will be visible 
from fhe Earfh, fhus causing a secular change in fhe observed pulse profile. This effecl has been ob¬ 
served in some fighl binary pulsars, for example, in PSRs B1913-I-16 (Weisberg el al. 1989; Kramer 
1998), B1534-I-12 (Slairs el al. 2004; Fonseca el al. 2014), J1141—6545 (Manchester el al. 2010), 
J1906-P0746 (Lorimer el al. 2006) and J0737-3039B (Perera el al. 2010). For PSRs B1534-P12 
and J0737—3039B we even have direcl measuremenfs of Ihe precession rates, which malch Ihe 
GR predictions (Slairs el al. 2004; Fonseca el al. 2014; Breton el al. 2008). The SKAl-MID will 
have comparable sensitivity to Ihe Arecibo telescope, meaning lhal signilicanl improvemenl in pre- 

^Effects due to the PN orders are those that are suppressed by the ratio of the orbital velocity, v, to the light speed, 
c, compared with those due to Newtonian gravity (Will 2014). It is denoted as “nPN” if the suppression factor is 
^(v2"/c2«). 
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Mass of Pulsar A (Mg^j^) 


Figure 2: The mass-mass diagram of PSR J0737—3039A/B, also known as the Double Pulsar (Kramer et 
ah, in prep.). In the figure the underlying gravitational theory is assumed to be GR. Shaded regions are 
forbidden by the individual mass functions because the sine of the orbital inclination must be < 1. The inset 
is an expanded view of the region of principal interest, where the intersection of seven curves at one point 
within measurement uncertainties proves the existence and uniqueness of a solution. For more details, see 
Kramer et al. (2006) and Breton et al. (2008). 


cession measurements of the pulsars that are in the Arecibo sky will await the SKA2. However, 
precessing pulsars not visible to the Arecibo telescope will be prime targets for SKAl-MID. For 
example, in the Double Pulsar, the precession rate of the B pulsar can be measured by modelling the 
orientation of its magnetosphere as it causes “flickering” eclipses of pulsar A (Breton et al. 2008). 
The rapidly changing eclipse-period flux density will be much better measured by both phases of 
the SKA, leading to a more stringent test of the precession rate’s agreement with GR. 

Since in GR the total angular momentum must be conserved up to 2PN order, a change in the 
direction of the pulsar spin must be compensated by a change in the direction of the orbital angular 
momentum, contributing a variation to the projected semimajor axis of the orbit (i / 0). This 
spin-orbit coupling is also known as the Lense-Thirring effect or frame-dragging effect, and it is a 
small effect that is numerically at 2PN order for slowly rotating neutron stars, whose measurement 
requires a high timing precision (Damour & Schaefer 1988). However, measuring it is important 
since it offers the potential to measure the moment of inertia of the pulsar (Damour & Schaefer 
1988; O’Connell 2004), of great interest for the EOS of cold dense nuclear matter (Watts et al. 
2014). 

For the Double Pulsar, the recycled pulsar (pulsar A) dominates the contribution to the Lense- 
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Thirring precession, because it is a faster rotator. The effect on i will be almost impossible to 
measure, not only because in this case there is a close alignment of the orbital angular momentum 
and the spin of pulsar A (Ferdman et al. 2013), but also because the system is viewed almost 
edge-on (Kramer et al. 2006). However, the SKA Galactic census will discover many NS-NS 
systems (scaling from the current population, we expect 100 from the SKAl and 180 from the 
SKA2 (Keane et al. 2014)), and there may well be binaries that will enable the direct measurement 
of X. In order to achieve the goal, a tight binary with good timing precision is needed. In addition, 
the misalignment angle between the pulsar spin and the orbital norm should be reasonably large, 
i.e., geodetic precession should be observable. This requires a relatively large “supernova kick 
velocity” for the second-formed compact object — in other words, an iron-core-collapse supernova, 
instead of an electron-capture supernova as likely happened in the Double Pulsar (Ferdman et al. 
2013). In such misaligned binaries, the frame dragging introduces nonlinear evolution in time in 
the longitude of periastron and its projected semimajor axis; thus one might be able to extract G) 
and X from radio timing. This would impose useful constraints on the system’s geometry (Wex & 
Kopeikin 1999). 

The Lense-Thirring effect also contributes to ft) with a magnitude comparable to that of ft)^™ 
(Damour & Schaefer 1988; Kramer & Wex 2009). For the Double Pulsar, this can in principle be 
extracted from the advance rate of periastron. To do this we must determine the component masses 
very precisely using two other precise PK parameters. One of them (the shape of Shapiro delay, 
s) already possesses the required precision (Kramer et al., in prep.). For the other one, Pj,, the 
orbital decay due to the emission of GW waves, we will need to estimate the distance to the pulsar 
carefully, in order to estimate precisely the kinematic corrections to it. This can only be done with 
the timing precision provided by the SKA. Again, if this Lense-Thirring contribution to ft) can be 
estimated precisely, we will then be able to estimate the moment of inertia of pulsar A with similar 
precision. Because this measurement will be done for a pulsar with an exquisitely well-determined 
mass, it will be a very important constraint on the EOS for dense matter (Lattimer & Schutz 2005; 
Kramer & Wex 2009; Watts et al. 2014). 

The SKA will also confirm general-relativistic orbital deformation for the first time. Damour 
& Deruelle (1986) used three different eccentricities to describe the timing of a pulsar. Two pa¬ 
rameters were introduced in the pulsar timing formula, namely 5e and 5r, where the former is in 
principle an observable (Damour & Taylor 1992). However, it only adds a tiny periodic variance 
to the TOAs and the signal has a strong correlation with the Einstein delay that is caused by the 
transverse Doppler effect and the gravitational redshift, its amplitude is js/X —e^ in binaries of 
mass ratio near unity, where e is the eccentricity of the orbit. Therefore, we will need highly ec¬ 
centric orbits to “amplify” the effect, and highly relativistic orbits with large periastron advance 
rate to break the degeneracy with the Einstein delay. The SKA is likely to find such a sysfem in fhe 
soufhern sky and fime if wifh high precision. We simulafe TOAs wifh orbifal paramefers similar fo 
fhaf of PSR B1913-I-16 (Hulse & Taylor 1975; Weisberg el al. 2010) wifh fhe TEMPO soflware^. If 
is assumed fhaf 60 TOAs wifh a precision ~ 1 /rs are obfained per monfh wifh fhe SKAl. After 10 
years of observation, fhe parameter can be defermined fo a precision ~ 10%. A more eccenfric 
(or a more relativistic) binary can achieve a beffer fesf. These resulfs supporf earlier simulations 


^http://tempo. sourceforge.net 
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done with the Double Pulsar (Kramer & Wex 2009). 

Where parameters in the timing model can be measured to high precision by other means, they 
can be restricted to an allowed range in the timing fit. This can lead to significantly reduced un¬ 
certainties for the remaining parameters in the presence of covariances. The most common source 
of external information is the estimation of position, proper motion and annual parallax using Very 
Long Baseline Interferometry (VLBI). The use of VLBI can improve gravitational tests when ap¬ 
plied to pulsar binaries. For instance, astrometric information provided by VLBI was used to cal¬ 
culate the kinematic contamination of the measured orbital period derivative of the Double Pulsar, 
improving the potential accuracy of tests of gravitational radiation to the 0.01% level (Deller et al. 
2009). In a second example, the previously unmeasured Shapiro delay for PSR 12222—0137 was 
precisely ascertained after the position, proper motion, and parallax were fixed using VLBI infor¬ 
mation (Deller et al. 2013). Since VLBI measures motions on the plane of the sky, it has different 
“blind spots” with respect to pulsar timing; VLBI has no “blind spots” to parallax at the ecliptic 
pole and proper motion in the ecliptic plane. With current VLBI instruments, precisions of 10 
micro-arcseconds for parallax and 10 micro-arcseconds per year for proper motion can be reached 
(Deller et al. 2013). The SKAl will allow a similar leap in VLBI precision as it will provide in the 
timing precision, and so VLBI using the SKA I will be able to deliver better results for parallax and 
proper motion than timing in most cases for nearby pulsars (Paragi et al. 2014). These measure¬ 
ments will be very important to subtract the Shklovskii contribution (which usually dominates the 
error budget), and the differential Galactic rotation, and vertical acceleration in the Galactic poten¬ 
tial (which also depend on relatively well-constrained Galactic models), thus measure the intrinsic 
orbital period change to high precision. 

3. Gravitational Dipole Radiation and Equivalence-Principle Tests 

One of the characteristics of GR is the property of “effacement” (Damour 1987). It states 
that the finite extents of two bodies in a binary are largely irrelevant to their orbital dynamics as 
long as they are well separated. Effacement is deeply rooted in the strong equivalence principle 
(SEP) (Will 1993). However, in alternative theories of gravity, the SEP is generally violated, and 
the strong gravitational fields associated with a NS may have measurable impact on the orbital 
dynamics of the binary and the trajectory of light propagation (Will 1993; Damour & Esposito- 
Earese 1992, 1993). Eigure 1 shows various binary systems, with the gravitational field strength 
characterised by the compactness and curvature of the orbit and the compactness of the companion 
star. 


3.1 Gravitational dipole radiation and time-variation of the gravitational constant 

A possible violation of the SEP could be identified by the occurrence of gravitational dipole ra¬ 
diation. The dipole radiation is at lower PN order than the canonical quadrupole radiation predicted 
by GR, thus in principle it can dominate over the quadrupole mode (Damour & Esposito-Earese 
1992, 1993; Will 1993, 2014). It carries away the gravitational energy of the orbit more rapidly 
than that in GR, collapsing the binary more quickly, and modifying the waveform of gravitational 
radiation that can be detected by future ground-based GW detectors. The leading change in the 
orbital period of a binary is given by the dipolar contribution with a form oc (ccj — 0 : 2 )^, 
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Figure 3: Constraints on the parameter space of scalar-tensor theories (figure courtesy of N. Wex). Coloured 
regions are excluded by current Solar System and pulsar experiments: Cassini spacecraft (Bertotti et al. 
2003), PSR J1738H-0333 (Freire et al. 2012b), and PSR J0348H-0432 (Antoniadis et al. 2013). The EOS is 
chosen as MPAl. Black dashed lines are based on simulations for a MSP-BH system with (Mbh,44msp) = 
(IOM 0 , 1 . 4 M 0 ), Ph=5d, and e = 0.8; the constraints from top to bottom are based on: 10 years with a 100- 
m class radio telescope, 5 years with the FAST telescope, and 5 years with the SKA (Liu et al. 2014). The 
vertical dashed purple line “JFBD” indicates the Jordan-Fierz-Brans-Dicke gravity. The horizontal dashed 
grey line “GAIA” indicates the limit expected from near-future Solar-System experiments, foremost from 
the astrometric satellite GAIA. 


where a,- (/ = 1,2) is the effective coupling strength between body i and the fields associated with 
the dipole moment. In the best-studied gravitational theory of this type, the scalar-tensor theory 
(Damour & Esposito-Farese 1992; Will 1993), we generally have Wns / otwD in a NS-WD binary 
(see the right panel of Figure 1). Therefore, NS-WD binaries are the ideal laboratories for testing 
these theories and, if they are correct, searching for the existence of gravitational dipole radiation, 
which if detected would falsify GR. 

In these NS-WD binaries, the combined radio observation of the pulsar and optical observation 
of the WD can give rise to precise measurements of two component masses (Antoniadis et al. 2012, 
2013; Wex 2014). The mass measurements can be used to predict the amount of gravitational 
radiation within a specific theory, and this in turn is confronted with the radio timing observation 
of the orbital decay caused by gravitational damping. As mentioned, the dipolar radiation can in 
principle be dominant over the quadrupolar radiation, therefore, agreement between the observed 
Pb and the quadrupole formula in GR stringently constrains the dipolar radiation contribution, as in 
the cases of PSRs J1738-I-0333 (Freire et al. 2012b) and J0348-I-0432 (Antoniadis et al. 2013). 

In Figure 3, the constraints from those two binaries are plotted in the tto-jSo parameter space 
mentioned in section 1. The best constraint from the Cassini spacecraft in the Solar System is also 
plotted. The limits from two binary pulsars stringently constrain the parameter space of scalar- 
tensor theories, especially when the jSq parameter is negative. The reason for this is that, as showed 
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by Damour & Esposito-Farese (1993), in this region NSs develop non-perturbative strong-field 
effeets like “spontaneous sealarisation” that are absent in the weak-field experiments in the Solar 
System; sueh effeets would eause a very large inerease in the emission of dipolar gravitational 
waves. The more massive the pulsar the less negative jSq has to be for the spontaneous sealarisation 
to appear. For this reason, the timing of PSR J0348-I-0432, a two-solar-mass NS in a relativistie 
orbit where the orbital deeay has been measured with some preeision, introduees uniquely stringent 
eonstraints on the allowed range of jSo (Antoniadis et al. 2013). 

In Jordan-Fierz-Brans-Dieke theory and (from a broader viewpoint) in the generie sealar- 
tensor theories, the sealar field plays fhe role of loeal gravifafional eonsfanf, implying fhaf fhe 
latter may vary wifh spaee and time (Damour & Esposito-Farese 1992; Will 1993). The possible 
time-variation in the gravitational eonstant eontributes to the deeay of a binary orbit (Damour et al. 
1988; Wex 2014). The extra eontribution ean be eonstrained by the measurement of Ph for binary 
pulsars by assuming a null eontribution from the dipole radiation. In general, however, a varying 
gravitational eonstant is aeeompanied by the dipole radiation (Fazaridis et al. 2009). With more 
than one binary pulsar, one ean also eonduet a joint eonstraint simultaneously on both the dipole 
radiation and the time-variation in the gravitational eonstant, utilising the faet that these two effeets 
have different dependenee on the orbital period of the binary (Fazaridis et al. 2009). The latest 
eonstraint from pulsars on the time-variation of the gravitational eonstant (Freire et al. 2012b), 
whieh relies on VFBl (Deller et al. 2008) and timing (Verbiest et al. 2008) of PSR J0437—4715 
and J1738-I-0333, is already eomparable to the best eonstraint from the Solar System experiments 
(Will 2014). Moreover, pulsar tests are sensitive to strong-field effeels on G (Wex 2014). 

These examples already show how imporfanf proper subfraefion of fhe eonfribufions from 
Galaefie aeeelerafion and fhe Shklovskii effeel is, and how erueial VFBl is for redueing uneerfain- 
fies in disfanees and proper motions. Furfhermore, in fhe near fulure a more aeeurafe undersfanding 
of fhe Galaefie pofenfial will also be available (e.g., by fhe GAIA safellife^), furlher improving fhe 
esfimafes of fhe kinemafie eorreefions of Pt,. 

The nexf generation of 30-m class opfical felescopes, for examples, fhe Thirfy Meier Telescope 
(TMT)^, fhe Gianf Magellan Telescope (GMT)^, and fhe European Exlremely Farge Telescope (E- 
EFT)’, may also significanlly advance fhe opfical observations of WDs, improving Iheir mass and 
dislance esfimafes; which will soon limil fhe precision of some of fhe currenf gravily fesls (Freire 
el al. 2012b). 

In shorl, in fhe SKA era, pulsar observations will improve significanlly and mainlain Iheir 
unique imporlance in lesling allernalive Iheories of gravily, bul will benefil exlensively from in¬ 
putting olher lechniques like VFBl and opfical observation. They will be crucial for fesls of some 
exfended gravifafional Iheories lhal incorporale dark matter and dark energy as gravifafional effeels. 

If any relalivislic PSR-BH binaries are discovered, Ihey will also be promising lesl beds for 
gravifafional dipolar radialion because fhe no-hair Iheorem also applies lo scalar-lensor Iheories of 
gravily; Ihis implies lhal BHs have no scalar charge^. A consequence of Ihis is lhal any radialive 

^http://sci.esa.int/gaia/ 

^http://www. tmt.org/ 

®http://www. gmto.org/ 

^http://www.eelt.org.uk/ 

^However, see Herdeiro & Radu (2014) for “Kerr BHs” with scalar hair after a complex, massive scalar field, 
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test that constrains (wns ~ c^bh) will therefore constrain Wns directly (Damour & Esposito-Farese 
1998). The constraints from simulations of a hypothetical PSR-BH system (Mbh = IOMq, Py = 
5 day and e = 0.8) are depicted as dashed lines in Figure 3 (Fiu et al. 2014). Five years of weekly 
observations with the SKA will constrain |o!o| down to the order of i^(10^^) (Fiu 2012; Wex 
et al. 2013), improving upon the expected constraint from the GAIA satellite. Worthy to note that, 
unlike NS-WD binaries, these tests with PSR-BH binaries will also constrain the parameter space 
with j3o ~ — 1. 

3.2 Equivalence-principle violation and its effects on orbital dynamics 

A possible violation of the SEP has also been investigated by tracing the orbital dynamics of 
small-eccentricity NS-WD binaries (Damour & Schaefer 1991). The differential free fall rates of 
the NS and the WD in the potential of the Milky Way will induce a characteristic time evolution 
of the orbital eccentricity vector (Damour & Schaefer 1991). The SEP violation tends to polarise 
the Faplace-Runge-Fenz vector of the orbit towards the direction of the Galactic acceleration. By 
using statistical studies on the eccentricity distribution of NS-WD binaries with a probabilistic 
assumption (Wex 2000), the dimensionless SEP violation parameter (the Nordtvedt parameter), T], 
was constrained to a comparable magnitude to that from the tests in the Solar System (Stairs et al. 
2005; Gonzalez et al. 2011; Will 2014). A direct test of SEP with binary pulsars that uses the 
time-variations in orbital parameters was also developed (Freire et al. 2012a). It does not rely on 
statistical or probabilistic assumptions. 

The recently discovered triple system, PSR J0337-I-1715 (Ransom et al. 2014), provides an 
opportunity to improve the SEP test by many orders of magnitude (Freire et al. 2012a). In this 
system, an outer WD orbits an inner NS-WD binary in less than one year. Usually, in the tests of 
SEP violation, the Galactic acceleration exerted on the binary is of the order i^(10^^*^ms^^), which 
is the force that would “polarise” the inner orbit in the case of SEP violation. In PSR J0337-I-1715, 
the outer WD provides a “polarising” acceleration that is more than one million times larger than 
the Galactic acceleration on the inner pair. As a rough estimate, we use the acceleration in Eq. (4) 
of Damour & Schaefer (1991), by replacing the relativistic rate of advance of periastron for the 
inner binary with the orbital angular velocity of the outer orbit (also see Eq. (8.18) in Will (1993)). 
Based on the current measurement precision of the orbital eccentricities in the triple system, a 
possible upper limit of the SEP violation parameter is estimated, | A| < 4 x 10^^, which is lO'^ times 
better than the current pulsar constraint (N. Wex, private communication). Due to this pulsar’s 
detectability with the Arecibo telescope, observations with SKAl-MID will mostly improve the 
SEP-violation test via the accumulation of more data^, but the SKA2 will improve the timing 
precision by a factor of ten, which should improve the test by a corresponding factor. 

It was estimated in Ransom et al. (2014) that there could be < 100 such systems with MSPs 
residing in the Galaxy. Since the SKA will discover almost all visible pulsars in the Galaxy (Keane 
et al. 2014), it is likely that it will discover other triple systems like PSR J0337-I-1715, and it is 
possible that the SKA might even discover a triple system consisting of an inner NS-WD close 

minimally coupled to GR, is included. 

®In principle, SKAl-MID is better due to less radio frequency interference (RFI), and much wider bands than that 
of the Arecibo telescope. 
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binary and where the outer star is another NS. Sueh a system eould be (again) many orders of 
magnitude more sensitive to SEP violation effeets than PSR J0337+1715. The reason is that the 
SEP violation effeets are proportional to (aj""®’' — (Ereire et al. 2012a), where 

are the eouplings of the inner bodies to the sealar field and is the eoupling of the outer body 
to the sealar field. If fhe oufer body is a WD, fhen ~ Oq (cCq is eoupling in fhe weak field), buf 
if if is a NS, fhen eould, aeeording fo some alfernafive fheories of gravify, be many orders of 
magnifude larger fhan ao- Eor sueh fheories, fhe prediefed SEP violafion effeefs would be roughly 
“otNs/oto” times larger than their predietion for PSR J0337+1715, therefore making the system 
mueh more sensitive to the SEP violation predieted in these theories. The non-deteetion of SEP 
violation in sueh a system would provide a stringent eonstraint on alternative theories of gravity 
that prediet sueh SEP violation. 

3.3 Local Lorentz invariance and local position invariance of gravity 

Besides the universality of the free fall, the other two aspeets of the Einstein equivalenee 
prineiple (Will 2014), namely loeal Eorentz invarianee (EEI) and loeal position invarianee (EPI) 
of gravity, ean also be tested with pulsar timing experiments. Seenarios with EEI violation have 
reeently raised great interests in the gravitational eommunity, for examples, in the Einstein-rether 
theory and the Hofava-Eifshitz theory. Two generie frameworks for testing deviations from GR 
are the parametrised post-Newtonian (PPN) formalism (Will 1993, 2014) and the standard-model 
extension (SME) (Bailey & Kosteleeky 2006), both of whieh eontain parameters for EEI violation. 

Generieally, violation of EEI and EPI leads to modifreations of the orbital dynamies of binary 
pulsars and the spin evolution of solitary pulsars in eharaeteristie ways (Damour & Esposito-Earese 
1992; Bailey & Kosteleeky 2006; Wex & Kramer 2007; Shao & Wex 2012; Shao et al. 2013; Shao 
2014; Shao & Wex 2013). The EEI-violating modifieations of the orbital dynamies of a binary 
pulsar introduee time-variations in the orbital eeeentrieity and orbital inelination, plus an extra 
eontribution to the periastron advanee rate, thus resulting in ehanges in the PK parameters, e, i, 
and d) (Damour & Esposito-Earese 1992; Bailey & Kosteleeky 2006; Shao & Wex 2012; Shao 
2014). The violation of EEI also leads to spin preeession of a solitary pulsar with respeet to a fixed 
direction^*’ (Shao ef al. 2013; Shao 2014), while fhe violafion of EPI leads fo spin preeession of a 
solifary pulsar around ifs “absolufe” aeeelerafion fowards fhe Galaefie eenfre (Shao & Wex 2013). 
They bofh ehange our view on fhe emission region of a pulsar, and resulf in eharaeferisfie ehanges 
versus lime in fhe pulse profiles. Inferesfingly, modifreafions of fhe orbilal dynamies and spin 
evolution are relaled and eonlrolled by one sel of paramelers in bofh PPN and SME frameworks. 
These modifieations ean be joinlly eonslrained and eross-eheeked. The eurrenl besl eonslrainls 
on fhe EEI of gravity and EPI of gravify are from pulsar experimenls (Sfairs el al. 2005; Shao & 
Wex 2012; Shao el al. 2013; Shao & Wex 2013; Shao 2014). If was shown fhal Ihese eonslrainls 
are proportional fo fhe timing preeision of binary pulsars. Wilh fhe SKA, fesfs will be improved 
dramalieally wilh heller sensilivilies. Simulations show fhal, if fhe signal fo noise ralio of pulse 
profiles improves by a faelor of 10 wilh fhe SKA, fhen fhe limils on Eorenlz-violaling paramelers 
improve fhe same faelor over fhe eurrenl besl ones wilh solely 10-year observalions (say, from 2020 

'^If there exists a preferred frame, the fixed direction will be the direction of the “absolute” velocity of the binary 
with respect to this preferred frame. 
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Figure 4: Characteristic pulsar timing residuals caused by the quadrupole moment of Sgr A* BH from 
simulated TOAs for three orbital phases (Liu et al. 2012) (figure courtesy of N. Wex). The orbital period is 
assumed to be Py = 0.3 yr and the eccentricity e = 0.5. The red line is the fitted model. 

to 2030). The acquisition of one stable profile biweekly is assumed for two MSPs in Shao et al. 
(2013). If these SKA observations are properly combined with the 20-year pre-SKA observations 
(say, from 2000 to 2020), the improvement will be a factor of ~ 50. In addition, new stable MSPs 
from the Galactic census will be very helpful in breaking the degeneracies of parameters in the test, 
due to the vectorial and/or tensorial field condensafions in fhe spacefime wifh such violafions (Shao 
2014). 

4. Pulsar - Black Hole Binaries 

BHs are predicted fo be fhe ulfimafe oufcome of massive sfars in fhe fheory of GR. Evidence 
for fhe existence of BHs in various asfrophysical environmenfs has accumulated during fhe pasf few 
decades, and is now rafher robusf (Narayan & McClinfock 2013). However, such observations are 
sfill indirecf, and no high-precision fesfs of BHs are yef possible. Tesfing BH physics is a fanlasfic 
and infriguing goal in fhe fields of asfronomy and physics. The goal of precision fesfs of properties 
of BHs will be achieved wifh fhe SKA (Kramer ef al. 2004). 

The combinafion of fhe sensifivify and fhe large field of view of fhe SKA will inevifably mul¬ 
tiply fhe number of pulsars by fens fo hundreds, making defecfion of rare objecfs like PSR-BH 
binaries possible. The shorfer infegrafion lime in pulsar search (say, a 10-min infegrafion as pro¬ 
posed) will be beneficial fo discoveries of exfremely relafivisfic binaries. A relafivisfic PSR-BH 
system in a clean environmenf represenfs a “holy grail” of pulsar asfronomy. A vasl unexplored pa¬ 
rameter space in ferms of gravifafional pofenfial, spacefime curvalure and objecl compacfness will 
Ihen become accessible (see Figure 1). Pulsar - sfellar-mass BH binaries, pulsar - intermediate- 
mass BH binaries, and pulsars orbifing fhe Sgr A* BH probe differenl new regimes of gravify where 
no ofher gravity experimenfs are yef able fo be performed wifh precision fesfs (Wex & Kopeikin 
1999; Pfahl & Loeb 2004). The pulsar - sfellar-mass BH binary may nafurally come from a pri- 
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mary binary composed of a > 10 Mq massive star and another > %Mq star, where the former ends 
up as a BH and the latter ends up as a NS. The pulsar - intermediate-mass BH binaries may exist 
in the centre of some globular clusters (Clausen et al. 2014; Hessels et al. 2014). Given the dense 
environment in the Galactic centre, there may be numerous NSs orbiting the Sgr A* BH, some of 
which may be detected as radio pulsars (Pfahl & Loeb 2004). However, no pulsars close enough 
to Sgr A* BH for gravity tests have been detected yet. Recently, a magnetar was detected through 
X-ray (Mori et al. 2013; Kennea et al. 2013) and radio (Eatough et al. 2013) observations, however, 
it is not close enough to the Sgr A* BH for precision tests of gravity (Eatough et al. 2013, 2014). 
Because of scattering by the interstellar medium, in general we will need the SKA2 to detect these 
pulsars in the Galactic centre (Wharton et al. 2012; Wex et al. 2013; Eatough et al. 2014), although 
this will depend on which bands are included in SKAl-MID, for example, the inclusion of Band 
5 will greatly enhance the prospects for doing this. Eor detecting the other two types of PSR-BH 
binaries, the SKAl will start to contribute significantly. 

In the case of relativistic PSR-BH binaries, the prominent effects from gravitation will, on one 
hand, reinforce better measurements of existing tests such as the Shapiro delay and the gravitational 
wave damping, and on the other hand, enable completely new tests of the “cosmic censorship 
conjecture” and the “no-hair theorem” (Damour & Esposito-Earese 1998; Wex & Kopeikin 1999; 
Eiu et al. 2012, 2014). The “cosmic censorship conjecture” states that for a realistic BH, there 
always exists an event horizon which prevents us from looking into its central singularity, thus 
preserving classical predictability in a spacetime. It imposes a maximum spin, S'max = GM^/c, on 
a rotating BH with a given mass M. The dimensionless spin of a BH, x = cSjGM^, should always 
satisfy the constraint % ^ ^ (Misner et al. 1973). In alternative gravitational theories, this bound 
may be slightly violated, as in the Einstein-Gauss-Bonnet theory, where, besides the Einstein- 
Hilbert term, the Gauss-Bonnet topological term is included in the Eagrangian of gravity (Kleihaus 
et al. 2011). The “no-hair theorem” states the uniqueness of a Kerr BH after it settles down. It 
relates the higher-order moments of a BH exclusively to its mass and spin. Eor example, if the 
“no-hair theorem” is valid, the dimensionless quadrupole of a BH, q = c^Q/G^M^ where Q is the 
quadrupole moment of BH, should satisfy q = —x^ (Misner et al. 1973). If a complex, massive 
scalar field, minimally coupled fo GR, is infroduced, fhen fhe “no-hair fheorem” could be violafed 
(Herdeiro & Radu 2014). 

Eor a relativistic PSR - Sgr A* BH system, fhe pulsar can be freafed as a fesf particle in fhe 
firsf-order approximafion (Hackmann & Eammerzahl 2008, 2012). The sfrafegy of performing fhe 
fesfs of fhe “cosmic censorship conjecfure” and fhe “no-hair fheorem” was invesfigafed in Eiu ef al. 
(2012). Eirsf of all, as usual in pulsar timing, we can estimate fhe mass of fhe Sgr A* BH fhrough 
fhe measuremenf of fhe periasfron advance rafe. After fhaf, a beffer mass determination and fhe 
geomefry information can be obfained wifh fhe defecfion of fhe Shapiro delay and fhe Einsfein 
delay. Wifh fhe help of fhe frame-dragging effecls (see secfion 2), fhe spin-orbif coupling fhen 
allows an inifial exfracfion of an upper limif of fhe spin magnifude of fhe Sgr A* BH. Thus an 
initial fesf of fhe “cosmic censorship conjecfure” becomes available by looking al fhe magnifude 
of X- This may enable us fo look info fhe bizarre possibilify of Sgr A* BH being a boson sfar 
(Kleihaus el al. 2012). Wifh a heller determination of fhe Sgr A* BH mass, fhe separation of fhe 
impacf on fhe periasfron advance from fhe BH spin is made plausible. This further delermines fhe 
3-dimensional spin of fhe BH wifh magnifude and direclion. An exlraclion of fhe periodic effecls 
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of the quadrupole from timing residuals will then allow a quantitative test of the “no-hair theorem”. 
The timing residual from the quadrupole of Sgr A* BH has a characteristic shape, as illustrated in 
Figure 4 for a hypothetical system with an orbital period Py ^ 0.3 yr and an eccentricity e ~ 0.5 
(Liu et al. 2012). According to simulations for an orbit with a period of several months, the spin 
and the quadrupole can be measured with precisions of ~ 0.1% and ~ 1% respectively after five 
years of observations with the SKA. Therefore, the test of the “no-hair theorem” can be performed 
at the ~ 1% level by then (Liu et al. 2012). 

Similarly, in the cases of stellar-mass BH and intermediate-mass BH binaries, the spin infor¬ 
mation of the BH can be extracted with the help of frame dragging as well if the orbit is tight 
enough (Wex & Kopeikin 1999). For the stellar-mass BH companion, a relativistic binary can pro¬ 
vide much better constraints on the scalar-tensor theories than the tests with gravitational waves 
from inspiralling binaries (Damour & Esposito-Farese 1998). These constraints will also be an 
order of magnitude better than the expected near-future constraint from the GAIA satellite when 
j3o is negative or small; see Figure 3 for results from a simulated PSR-BH system (Liu et al. 2014). 

5. Conclusion and Outlook 

By the 2020s, radio astronomy will undergo a revolution brought by the onset of the SKAl. Its 
pulsar survey will deliver a large population of pulsars, among which there will be ~ 1800 MSPs. 
The most stable ones of them will be used to perform a number of precision gravity tests. The 
improvement in the sensitivity of the timing precision allows us a closer look at many aspects of 
the foundation of gravity. In the era of the SKA, the existing tests will be pushed to new extremes, 
and completely new tests will also be made possible with systems including NS-NS binaries, NS- 
WD binaries and NS-BH binaries. Advances from the theoretical side will also contribute to the 
area of experimental gravity with novel tests. It is indeed a virtue of nature that pulsar binaries 
of different types are equipped with different characteristics to perform different tests with. For 
example, the NS-NS systems provide extremely precise tests of GR, in particular periodic, non- 
dissipative effects in orbital dynamics {e.g. with the Einstein delay and general-relativistic orbital 
deformation), while the NS-WD and NS-BH systems are suitable to study the dissipative radiative 
property of gravity in the presence of asymmetry in gravitational binding energy {e.g. tests of the 
gravitational dipole radiation and the constancy of the gravitational constant). Eurthermore, PSR- 
BH binaries will be useful to test the BH physics {e.g. the “cosmic censorship conjecture” and the 
“no-hair” theorem). Tests of gravitational theories from many facets in the era of the SKA will 
profoundly advance our understandings of fundamental physics. 

The output of gravity tests with pulsars has two possibilities. In the first possibility, GR is 
falsified. This will certainly provide an important hint for our further understanding of the funda¬ 
mental interaction, gravity, and the nature of spacetime. In the second possibility, GR is confirmed 
to new precision with flying colours. This will further reinforce our confidence in using GR as a 
tool to investigate other phenomena, for example, using GR effects to make mass measurements 
of binary pulsars, helping to deepen our understandings of non-perturbative QCD effects (Watts 
et al. 2014). It also casts experimental constraints on constructing new gravitational theories in 
relevant areas, for example, to explain the nature of dark matter and dark energy. In both cases. 
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tests of gravitational theories with great precision will provide important clues to reconcile GR and 
quantum principles. 

As a final remark, with the better timing capability and new discoveries of stable pulsars from 
the SKA, the pulsar timing array (PTA) projects will leap towards a new era as well. Comple¬ 
mentary to the strong-held tests discussed here, PTAs make tests of gravitational theories with 
gravitational waves in the radiative regime available (Janssen et al. 2014). 
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